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ABSTRACT 

Aims. We present a multi-frequency study of the supernova remnant MCSNR 10512-6707 in the Large Magellanic Cloud. 

Methods. We used new data from XMM-Newton to characterise the X-ray emission and data from the Australian Telescope Compact 
Array, the Magellanic Cloud Emission Line Survey, and Spitzer to gain a picture of the environment into which the remnant is 
expanding. We performed a morphological study, determined radio polarisation and magnetic field orientation, and performed an 
X-ray spectral analysis. 

Results. We estimated the remnant’s size to be 24.9 (±1.5) x 21.9 (±1.5) pc, with the major axis rotated ~ 29° east of north. Radio 
polarisation images at 3 cm and 6 cm indicate a higher degree of polarisation in the northwest and southeast tangentially oriented to 
the SNR shock front, indicative of an SNR compressing the magnetic field threading the interstellar medium. The X-ray spectrum 
is unusual as it requires a soft (~ 0.2 keV) collisional ionisation equilibrium thermal plasma of interstellar medium abundance, in 
addition to a harder component. Using our fit results and the Sedov dynamical model, we showed that the thermal emission is not 
consistent with a Sedov remnant. We suggested that the thermal X-rays can be explained by MCSNR J0512-6707 having initially 
evolved into a wind-blown cavity and is now interacting with the surrounding dense shell. The origin of the hard component remains 
unclear. We could not determine the supernova type from the X-ray spectrum. Indirect evidence for the type is found in the study of 
the local stellar population and star formation history in the literature, which suggests a core-collapse origin. 

Conclusions. MCSNR J0512-6707 likely resulted from the core-collapse of high mass progenitor which carved a low density cavity 
into its surrounding medium, with the soft X-rays resulting from the impact of the blast wave with the surrounding shell. The unusual 
hard X-ray component requires deeper and higher spatial resolution radio and X-ray observations to confirm its origin. 
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1. Introduction 

Supernova remnants (SNRs) are formed via the interaction of 
ejecta from violent supernova (SN) explosions with the sur¬ 
rounding medium. The supernovae (SNe) can result from either 
the core-collapse (CC) of a massive star or from the explosion 
of a carbon-oxygen white dwarf in a binary system which ex¬ 
ceeds the Chandrasekhar limit, either via accretion from a stel¬ 
lar donor or white dwarf mergers, so-called Type la explosions. 
The SN ejecta blasted into the interstellar medium (ISM) con¬ 
tain heavy elements produced in the stellar progenitors and via 
explosive nucleosynthesis which chemically enriches the ISM. 
The energy deposited by the explosions drive the mechanical 
evolution of the ISM and cosmic rays can be accelerated in the 
fast moving shocks (see Vink 2012, for a review). Studies of 
these objects in the Milky Way are problematic due to uncertain 


* Based on observations obtained with XMM-Newton, an ESA sci¬ 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA 


distance estimates and high foreground absorption, which par¬ 
ticularly affects evolved SNRs whose low temperature plasmas 
are extremely susceptible to foreground absorption. 

The Large Magellanic Cloud (EMC) is an excellent galaxy 
for the study of SNRs due to its ideal observational properties. 
At a distance of 50 kpc (di Benedetto 2008) it is sufficiently 
close that its stellar population and diffuse structure is resolved 
in most wavelength regimes. The EMC is almost face-on (incli¬ 
nation angle of 30^0°, van der Marel & Cioni 2001; Nikolaev 
et al. 2004) and the modest extinction in the line of sight (av¬ 
erage Galactic foreground Ah ~ 7 x 10^® cm“^) means optical 
and X-ray observations of SNRs are only slightly affected by 
foreground absorption, whereas its location in one of the coldest 
parts of the radio sky (Haynes et al. 1991) allows for improved 
radio observations without interference from Galactic emission. 
The total number of confirmed SNRs in the EMC now stands at 
~ 60 (Maggi et al. 2015; Bozzetto et al. 2015, and references 
therein). Recent additions include Bozzetto et al. (2014), Maggi 
et al. (2014), Kavanagh et al. (2015a), Kavanagh et al. (2015b). 
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Typically, objects are classified as SNRs based on satisfying 
certain observational criteria. The Magellanic Cloud Supernova 
Remnant (MCSNR) Database^ requires that at least two of the 
following three observational criteria must be met: significant 
Ha, [S It], and/or [O iii] line emission with an [S ii]/Ha flux ratio 
> 0.4 (Mathewson & Clarke 1973; Fesen et al. 1985); extended 
non-thermal radio emission; and extended thermal X-ray emis¬ 
sion. A discussion of the significance of each of these classifica¬ 
tion criteria is given in Filipovic et al. (1998). 

The ROSAT catalogue of X-ray sources in the LMC of Haberl 
& Pietsch (1999) contained the source [HP99] 483 at a J2000 
position of RA = 05‘'12“28.0'* and Dec = -61^01'^2T which is 
located in the wider DEM L97 HII region (Davies et al. 1976). 
The source is also projected near the edge of the molecular cloud 
[WHO2011] A126 (Wong et al. 2011) to its south and southwest. 
However, due to the faintness of the source, only coarse proper¬ 
ties of [HP99] 483 were determined by Haberl & Pietsch (1999) 
and the source did not make it into the final list of ROSAT de¬ 
tected SNRs and candidate SNRs. 

The study of Reid et al. (2015) identified the optical coun¬ 
terpart to [HP99] 483 using deep Ha imaging of the LMC per¬ 
formed with the United Kingdom Schmidt Telescope (UKST). 
Spectroscopic follow-up of several regions in the object with 
the Anglo Australian Telescope and the 1.9 m telescope at the 
South African Astronomical Observatory revealed the charac¬ 
teristically high [S ii]/Ha values, consistent with an SNR. The 
radio counterpart was also detected using data from the Aus¬ 
tralia Telescope Compact Array (ATCA), whose spectral index 
of -0.52 (+0.04) is consistent with with the typical value of -0.5 
for non-thermal emission from SNRs. Therefore, the object was 
confirmed as an SNR by Reid et al. (2015), who assigned the 
identifier MCSNR J0512-6707. A study of the surrounding dust 
emission using Spitzer data was also presented, with the authors 
noting that the SNR is located in a complex dust environment 
with the SNR appearing to be encased in dust, notably by a bright 
filament in the east and north. Using optical imaging, these au¬ 
thors determined the size of the remnant to be 13.5 x 15.5 pc and 
estimated the age to be 2-5 kyr assuming that the remnant is in 
the Sedov phase. 

In this paper we present new X-ray observations of MC¬ 
SNR J0512-6707 with XMM-Newton. As we will show, the new 
X-ray data revealed that the size of MCSNR J0512-6707 is sig¬ 
nificantly larger (24.9 (±1.5) x 21.9 (+1.5) pc, see Section 4.1) 
than determined by Reid et al. (2015), with the optical emission 
only present in the northeast of the remnant. In addition, we will 
show that MCSNR J0512-6707 is likely not in the Sedov phase 
(Section 4.2). We performed a similar study to assess the mor¬ 
phology and physical characteristics in light of these new find¬ 
ings, making use of the data from the Magellanic Cloud Emis¬ 
sion Line Survey (MCELS, Smith et al. 2006), radio data from 
ATCA, and IR data from Spitzer reported by Reid et al. (2015). 

Our work on the multi-frequency study of MC¬ 
SNR J0512-6707 using X-ray, radio, IR, and optical emission 
line data is arranged as follows: the observations and data 
reduction are described in Section 2; data analysis is outlined 
in Section 3; results are given and discussed in Section 4; and 
finally we summarise our work in Section 5. 


' http://www.mcsnr.org/ 


2. Observations and data reduction 

2.1. X-ray 

XMM-Newton (Jansen et al. 2001) observed MC¬ 
SNR J0512-6707 on May 17 2014 (Obs. ID 0741800201, 
PI P. Kavanagh). The primary instrument for the observation 
was the European Photon Imaging Camera (EPIC), which 
consists of a pn CCD (Striider et al. 2001) and two MOS CCD 
(Turner et al. 2001) imaging spectrometers. However, a “false 
radiation” warning due to a flaring optical source in the EOV 
affected the EPIC-MOS CCDs during this observation, resulting 
in both EPIC-MOS cameras being shut for extended periods. 
Because of this, the observation was repeated and a second 
dataset was obtained on June 8 2014 (Obs. ID 0741800301, PI 
P. Kavanagh). We reduced all available data (i.e., the usable data 
from Obs. ID 0741800201 and Obs. ID 0741800301) using the 
standard reduction tasks of SAS^ version 14.0.0, filtering for 
periods of high particle background. This resulted in combined 
flare-filtered exposure times of ~ 42 ks for EPIC-pn, ~ 48 ks for 
EPIC-MOS 1, and ~ 46 ks for EPIC-MOS2. 

2.2. Optical 

The MCELS observations (Smith et al. 2006) were taken with 
the 0.6 m University of Michigan/Cerro Tololo Inter-American 
Observatory (CTIO) Curtis Schmidt Telescope equipped with 
a SITE 2048 x 2048 CCD, producing individual images of 
1.35° X 1.35° at a scale of 2.3" pixeU'. The survey mapped 
both the LMC (8° x 8°) and the SMC (3.5° x 4.5°) in narrow 
bands covering [Oiii]/I5007 A, Ha, and [Sii]/I6716, 6731 A, in 
addition to matched green and red continuum bands. The survey 
data were flux calibrated and combined to produce mosaicked 
images. We extracted cutouts centred on MCSNR J0512-6707 
from the MCELS mosaics. We subtracted the continuum images 
from the corresponding emission line images, thereby remov¬ 
ing the stellar continuum and revealing the full extent of the 
faint diffuse emission. Einally, we divided the continuum sub¬ 
tracted [S It] image by the continuum subtracted Ha image to 
get a [Sii]/Ha map of MCSNR J0512-6707, with regions of 
[S ii]/Ha > 0.4 indicative of the presence of an SNR (Mathewson 
& Clarke 1973; Eesen et al. 1985). 

2.3. Radio 

2.3.1. Radio continuum 

MCSNR J0512-6707 was observed with the Australian Tele¬ 
scope Compact Array (ATCA) on November 15 and 16 2011 
(project C634), using the new Compact Array Broadband Back¬ 
end (CABB). The ATCA array configuration EW367 was used 
and observations were taken simultaneously at T = 3 and 6 cm 
(v-9 000 and 5 500 MHz) using the dual-frequency mode. Base¬ 
lines formed with the 6* ATCA antenna were excluded as the 
other five antennas were arranged in a compact configuration. 
The observations were carried out in the so called “snap-shot” 
mode, totalling ~50 minutes of integration over a 14 hour pe¬ 
riod. PKS B1934-638 was used for flux density calibration and 
PKS B0530-727 was used for secondary (phase) calibration. The 
phase calibrator was observed twice every hour for a total 78 
minutes over the whole observing session. The miriad^ (Sault 


^ Science Analysis Software, see http://xmiti. esac. esa. int/sas/ 
^ http://WWW.atnf.csiro.au/computing/software/miriad/ 
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et al. 1995) and karma"^ (Gooch 1995) software packages were 
used for reduction and analysis. 

The relatively small size of the SNR coupled with the observ¬ 
ing array focusing on the shorter baselines, led to the remnant re¬ 
maining unresolved from these observations, and therefore, the 
radio morphology of the remnant cannot be readily seen. As an 
alternative, we use archival observations from the ATCA (project 
C1395) at d = 20 cm (v = 1384 MHz; bandwidth = 128 MHz). 
These observations were taken using ATCA arrays 1.5A and 6B 
on April 18 2005 and June 24 2005, respectively, for a com¬ 
bined integration time of 819.7 minutes. As these arrays focus 
on longer baselines, we were able to resolve the SNR, as seen in 
Fig. 1. 


2.3.2. CO 

To gain an understanding of the molecular environment in the re¬ 
gion of MCSNR J0512-6707, we used data from the Magellanic 
Mopra Assessment (MAGMA) survey by Wong et al. (2011)^. 
MAGMA is the mapping survey of the Large and Small Mag¬ 
ellanic Clouds in the CO(l-O) line using the 22 m Mopra Tele¬ 
scope of the Australia Telescope National Facility. The survey 
data have a spatial resolution of ~ 45" corresponding to 11 pc at 
the LMC distance. 
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Fig. 1. 20 cm radio-continuum image of MCSNR J0512-6707 overlaid 
with contours at 1.5, 2.0, 2.5, 3.0, and 3.5 mjy/beam. The ellipse in the 
lower left comer represents the synthesised beamwidth of 13.0"x 12.1". 


2.4. Infrared 

The cold environment surrounding MCSNR J0512-6707 can be 
revealed by infrared (IR) emission. To aid in the discussion of 
the morphology and environment of the remnant, we make use 
of data from the SAGE survey of the LMC (Meixner et al. 2006) 
with the Spitzer Space Telescope (Werner et al. 2004). During the 
SAGE survey, a 7° x 7° area of the LMC was observed with the 
Infrared Array Camera (IRAC, Eazio et al. 2004) in the 3.6 pm, 
4.5 pm, 5.8 pm, and 8 pm bands, and with the Multiband Imag¬ 
ing Photometer (MIPS, Rieke et al. 2004) in the 24 pm, 70 pm, 
and 160 pm bands. The MIPS 24 pm images provide us with 
a picture of the stochastically, thermally, and radiatively heated 
dust in the region of MCSNR J0512-6707, with spatial res¬ 
olution comparable to XMM-Newton, to give an indication of 
the distribution of cool material. We obtained the 24 pm MIPS 
mosaicked, flux-calibrated (in units of MJy sr ') images pro¬ 
cessed by the SAGE team from the NASA/IPAC Infrared Sci¬ 
ence Archive®. The pixel sizes correspond to 4.8" for the 24 pm 
band. 


3. Analysis 

The object was detected in radio, optical, and X-rays as an ex¬ 
tended source with a size of ~L7 x L5 (see Section 4.1). Here 
we describe the analysis of the SNR in these multi-wavelength 
datasets. 

3.1. Radio 

Linear polarisation images of MCSNR J0512-6707 were 
created at the same resolution at 6 and 3 cm using the Q and 
U Stokes parameters from the miriad task IMPOL. The mean 
fractional polarisation was calculated using flux density and 

http: //WWW. atnf.csiro.au/computing/software/karma/ 

^ See also http://inmwave. astro. Illinois. edu/magma/ 

® See http : //irsa.ipac.caltech.edu/data/ SPITZER/SAGE/ 


polarisation: 


- ( 1 ) 

O/ 

where Sq,S u, and 5/ are integrated intensities for the Q, U, and 
I Stokes parameters. Mean polarisation across the remnant was 
found to be 12.0 (+1.4)% at 6 cm and 7.1 (+1.9)% at 3 cm. Eig. 2 
shows the polarisation for both wavelengths, where a signal-to- 
noise cut-off of 2cr was used for the Q and U images, and a level 
of 6cr for the intensity image. Values that fell below these cut-off 
levels are blanked in the output image. As we detected polarised 
emission at two frequencies, we attempted to determine the Eara- 
day rotation through the miriad task imrm. However, as there is 
not any notable rotation between the vectors in the two images, 
the resulting errors were in the same order as the detected rota¬ 
tion measure. 

3.2. X-ray imaging 

We produced images and exposure maps in various energy bands 
from the flare-filtered event lists for each EPIC instrument in 
each observation. We filtered for single and double-pixel events 
(PATTERN < 4) from the EPIC-pn detector, with only single pixel 
events considered below 0.5 keV to avoid the higher detector 
noise contribution from the double-pixel events at these ener¬ 
gies. All single to quadruple-pixel events (PATTERN < 12) were 
considered for the MOS detectors. 

We used three energy bands suited to the analysis of the spec¬ 
tra of SNRs. A soft band from 0.3-0.7 keV includes strong lines 
from O; a medium band from 0.7-1.1 keV comprises Ee L-shell 
lines as well as Ne Hea and Lya lines; and a hard band (1.1- 
4.2 keV) which includes lines from Mg, Si, S, Ca, Ar, and non- 
thermal continuum if present. 

We subtracted the detector background from the images us¬ 
ing filter-wheel-closed data (EWC). The contribution of the de¬ 
tector background to each EPIC detector was estimated from the 
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Fig. 2. Polarisation of MCSNR J0512-6707 at 6 cm (top) and 3 cm 
(bottom) overlaid on an intensity image at the same wavelength. The 
6 cm contours are 1 to 11 mjy in steps of 2 mjy, and the 3 cm contours 
are 0.9 to 9.9 mJy in steps of 1.8 mJy. The ellipse in the lower left 
corners represents the synthesised beamwidth of 37.5"x 25.0", and the 
line directly below the ellipse represents a polarisation vector of 100 per 
cent. 


count rates in the corner of the images, which were not exposed 
to the sky. We then subtracted appropriately-scaled FWC data 
from the raw images. We merged the EPIC-pn and EPIC-MOS 
images from each observation into combined EPIC images and 
performed adaptive smoothing of each using an adaptive tem¬ 
plate determined from the combined energy band (0.3^.2 keV) 
EPIC image. The sizes of Gaussian kernels were computed at 
each position in order to reach a signal-to-noise ratio of five, 


setting the minimum full width at half maximum of the kernels 
to ~ 5". In the end the smoothed images were divided by the 
corresponding vignetted exposure maps. Einally, we produced a 
three-colour image of MCSNR J0512-6707 which is shown in 
Pig. 3-top left. 

3.3. X-ray spectral analysis 

Por the spectral analysis, we made use of the EPIC-pn and EPIC- 
MOS data. However, because of a low number of total counts, 
the MOS data from Ohs ID 0741800201 were omitted from 
our analysis. We extracted source and background spectra from 
vignetting-weighted event lists for each EPIC instrument which 
include a correction for the effective area variation across the 
source and background. This was achieved using the SAS task 
evigweight. The SNR spectra were extracted from elliptical 
regions encompassing the X-ray extent of the SNR. To ensure 
enough counts for a solid characterisation of its components, the 
background spectra were extracted from a larger annulus sur¬ 
rounding the SNR with point sources excluded. All spectra were 
rebinned so that each bin contained a minimum of 30 counts 
to allow the use of the statistic during spectral fitting. The 
EPIC-pn and EPIC-MOS source and background spectra from 
each observation were fitted simultaneously using XSPEC (Ar- 
naud 1996) version 12.8.2p with abundance tables set to those 
of Wilms et al. (2000), photoelectric absorption cross-sections 
set to those of Balucinska-Church & McCammon (1992), and 
atomic data from ATOMDB 3.0.1^ with the latest equilibrium 
and non-equilibrium data. 

3.3.1. X-ray background 

Detailed descriptions of the X-ray background constituents and 
spectral modelling can be found in Bozzetto et al. (2014) and 
Maggi et al. (2014). Here we briefly summarise the treatment of 
the X-ray background in the case of MCSNR JOS 12-6707. 

The X-ray background consists of the astrophysical X-ray 
background (AXB) and particle induced background. The AXB 
typically comprises four or fewer components (Snowden et al. 
2008; Kuntz & Snowden 2010), namely the unabsorbed ther¬ 
mal emission from the Local Hot Bubble, absorbed cool and hot 
thermal emission from the Galactic halo, and an absorbed power 
law representing unresolved background active galactic nuclei 
(AGN). The spectral properties of the background AGN com¬ 
ponent were fixed to the well known values of T ~ 1.46 and a 
normalisation equivalent to 10.5 photons cm“^ s ' sr"* * at 1 keV 
(Chen et al. 1997). The foreground absorbing material comprises 
both Galactic and EMC components. The foreground Galactic 
absorption component was fixed at 5.1 x 10^" cm“^ based on 
the Dickey & Lockman (1990) HI maps, determined using the 
HEASARC Ah Tool^, while the foreground EMC absorption 
component, with abundances set to those of the EMC, was al¬ 
lowed to vary in the spectral fits. 

The particle-induced background of the EPIC consists of the 
quiescent particle background (QPB), instrumental fluorescence 
lines, electronic read-out noise, and residual soft proton (SP) 
contamination. To determine the contribution of these compo¬ 
nents we made use of vignetting corrected EWC data. We ex¬ 
tracted EWC spectra from the same detector regions as the ob¬ 
servational source and background spectra. The EPIC-pn and 

’ http://www.atomdb.org/index.php 

* http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/ 

w3nh.pl 
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Fig. 3. Top left: XMM-Newton EPIC image of MCSNR JOS 12-6707 in false colour with RGB corresponding to 0.3-0.7 keV, 0.7-1. 1 keV, and 
1.1-4.2 keV. The image is overlaid with [Sii]/Ha contours with the lowest level corresponding to [Siil/Hor = 0.4, and the remaining levels at 
25%, 50%, and 75% of the maximum above this level. Top right: Continuum subtracted MCELS image of MCSNR 10512-6707 with Ho- in 
red and [S ii] shown in green overlaid with 0.3-0.7 keV contours. The lowest contour level represents 3cr above the average background surface 
brightness, with the remaining levels marking 25%, 50%, and 75% of the maximum above this level. Bottom left: Same as top left but with 20 cm 
radio contours from Eig. 1 overlaid. The contours are at 1.5, 2.0, 2.5, 3.0, and 3.5 mjy/beam. Bottom right: Spitzer MIPS 24 image of the 
MCSNR J0512-6707 region with the X-ray contours from top right in white and [Sn]/Htt contours from top left in green. The cyan contours 
represent the CO emission from the MAGMA survey and delineate the location of the molecular cloud [WHO2011] A126. The lowest level 
corresponds to 1.2 K km s“* (the approximate sensitivity limit of the survey), increasing in steps of 1 K km s“*. The image scale is the same as in 
all other panels. 


EPIC-MOS FWC spectra were fitted with the empirical mod¬ 
els developed by Sturm (2012) and Maggi et al. (2015), respec¬ 
tively. Since these spectral components are not subject to the in¬ 
strumental response, we used a diagonal response in XSPEC. 


The resulting best-fit model was included and frozen in the fits 
to the observational spectra, with only the widths and normalisa¬ 
tions of the fluorescence lines allowed to vary. We also included 
a multiplicative constant to normalise the continuum to the ob- 
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servational spectra using the high energy tail (E > 5 keV) where 
the QPB component dominates. The residual SP contamination 
was fitted by a power law not convolved with the instrumental 
response (Kuntz & Snowden 2008), which was only required in 
the spectra of Obs ID 0741800301. 

3.3.2. Source emission 

Depending on the age and ambient ISM density of an SNR, the 
X-ray emission is expected to be dominated by shock-heated 
ISM swept-up by the blast wave and/or reverse shock-heated 
ejecta. The swept-up ISM component should have abundances 
consistent with the LMC, whereas the ejecta component will be 
characterised by strong emission lines from its metal content. 
A cursory scan of the extracted spectra showed emission lines 
below 1 keV, however, these did not appear to be significantly 
enhanced and do not suggest a shocked-ejecta component. The 
spectrum above 1 keV was notably featureless with no evidence 
for emission lines. 

As a first attempt to model the X-ray emission from MC- 
SNR 10512-6707, we included a single thermal plasma model 
of LMC abundance absorbed by foreground Galactic and LMC 
material on top of the AXB and particle-induced background. 
We initially assumed that the plasma was in collisional ionisa¬ 
tion equilibrium (CIE), and therefore applied the vapec (Smith 
et al. 2001) model. The abundance of the LMC absorption com¬ 
ponent was fixed to the LMC values (0.5 Zq, Russell & Dopita 
1992). The best-fit vapec model (see Table 1 for fit results), 
with kT ~ 0.2 keV, adequately accounts for the emission below 
~ 1 keV but cannot explain the harder emission above 1 keV 
(see Fig. 4-left). This is reflected in the fit-statistic with reduced 
(xl) - 1-31- In initial test fits we also allowed the O and Fe 
abundance parameters in the vapec component to vary to search 
for signatures of ejecta emission, with O and Fe being the dom¬ 
inant constituents of CC and Type la explosive nucleosynthesis 
products, respectively. However, the resulting confidence ranges 
on these parameters were consistent with the LMC values, and 
we therefore kept the fixed LMC abundances. 

We relaxed our assumption of a CIE plasma and made use 
of a non-equilibrium ionisation (NEI) model in XSPEC ap¬ 
propriate for SNRs, namely the plane-parallel vpshock model 
(Borkowski et al. 2001). The vpshock model features a lin¬ 
ear distribution of ionisation ages behind the shock which is 
more realistic than single ionisation age models such as vnei. 
While this model resulted in a substantially improved fit-statistic 
(xl = 1.11, the best-fit model suggested a very high plasma 
temperature {kT > 4 keV) and very low ionisation parameter 
(t„ < 10'® s cm“^). Such spectral parameters are necessary to 
account for the hard tail above 1 keV. Allowing abundance pa¬ 
rameters to vary to account for possible ejecta signatures did not 
improve the fit. Given the very high plasma temperature and very 
low ionisation parameter, we were reluctant to take the fit results 
as the true physical description of the SNR emission, though it 
was clear that some component must be introduced to account 
for the harder emission above 1 keV. 

Therefore, we applied a simple two component fit with a 
thermal plasma and an additional power law (powerlaw in 
XSPEC). In trial fits we applied both the vapec and vpshock 
model to represent the softer emission component. However, 
since the plasma temperature and ionisation parameter (though 
poorly constrained) are consistent with a plasma in CIE, we pro¬ 
ceed using the vapec model to account for the soft component. 
The vapec+powerlaw model resulted in an improved fit over 
both single thermal plasma models WXhxl - L04. The vapec 


model parameters are largely consistent with the single vapec 
fit {kT ~ 0.2 keV), and the harder emission is accounted for 
by the power law with E = 2.07 (1.81 - 2.26), see Table 1 and 
Fig. 4-right. We note here that we also tried a two-temperature fit 
employing the vpshock model for the harder component. While 
this yielded a good fit {Xv - LIO) the best-fit temperature and 
ionisation parameter of the vpshock component were poorly 
constrained with the best-fit driven to high temperature and low 
ionisation parameter to handle the hard, featureless emission. 
However, we note that the hard emission may be modelled with 
either a power law or a thermal plasma. 

4. Results 

4.1. Multi-wavelength morphology 

The X-ray emission from MCSNR J0512-6707 is approx¬ 
imately circular, though slightly extended in the northeast- 
southwest direction, with the brightest region of emission in the 
northeast (Fig. 3 top-left). Emission from the remnant is detected 
in all energy bands with the softest emission tracing the outer 
edge, and the hardest emission localised in the northeastern re¬ 
gion and correlated with the brightest radio emission (Figs. 3 
bottom-left and 1) and the enhanced [Sii] emission detected in 
the optical (Figs. 3 top-left and top-right). The MIPS 24 pm 
image (Fig. 3 bottom-right) reveals a possible reason for this. 
The remnant is projected against a denser eastern region, as 
evident by a large arm of cool material emanating from the 
[WHO2011] A126 molecular cloud to the south and southwest, 
bordering the eastern front of the SNR, as noted by Reid et al. 
(2015). If indeed the SNR is expanding into this cool material to 
the east, the higher density in the east compared to the west can 
explain the observed multi-wavelength morphological properties 
of MCSNR J0512-6707. There is no evidence that the SNR is 
interacting with [WHO2011] A126 itself. 

The 20 cm radio morphology exhibits a more shell-like struc¬ 
ture than the X-ray or optical. The radio polarisation images 
(Fig. 2), the 6 cm in particular, indicate a higher degree of po¬ 
larisation in the northwest and southeast, tangentially oriented to 
the SNR shock front. Such a polarisation structure typically in¬ 
dicates evolved SNR shells (Fiirst & Reich 2004), with the tan¬ 
gentially oriented magnetic fields indicative of compression of 
the ISM threading magnetic field. 

To estimate the size of the X-ray remnant we determined 
the average background surface brightness and corresponding 
standard deviation (cr) in the 0.3-0.7 keV band. We defined the 
edge of the SNR as regions where the extended emission sur¬ 
face brightness rises to 3cr above the average background and 
fitted an ellipse to this contour. The error on the fit was deter¬ 
mined by quantifying the standard deviation of points on the 
contour from the best-fit ellipse. We determined a best-fit el¬ 
lipse centred on the J2000 coordinates RA = 05*’12™26.87'* and 
Dec = -67‘'07”2L7^ of size F72 (+0a0) x L51 (+0a0), cor¬ 
responding to 24.9 (+1.5) X 21.9 (±1.5) pc at the LMC distance, 
with the major axis rotated ~ 29° east of north. The best-fit di¬ 
mensions and error are shown in Fig. 5. 

4.2. X-ray emission 

The X-ray spectrum of MCSNR J0512-6707 is rather puzzling. 
There are obvious thermal features below ~ 1 keV, which are 
easily fitted with a thermal plasma component in CIE. However, 
this leaves significant residuals above ~ 1 keV. An NEI plasma 
can account for both the soft and hard parts of the spectrum, but 
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Fig. 4. XMM-Newton EPIC-pn spectrum of MCSNR J0512-6707 from Obs. ID 0741800201, shown alone for clarity. In each panel the best-fit 
model is given by the solid black line, the purple dash-dot lines mark the AXB components, the blue dash-dot-dot-dot line shows the combined 
contributions of the QPB, residual SPs, instrumental fluorescence lines, and electronic noise. Left: Best-fit vapec model shown by the dashed 
red line. Right: Best-fit vapec+powerlaw model shown by the dashed red line and dashed green line, respectively. The model fit parameters are 
given in Table 1. Note that the inshumental fluorescence line normalisations are free in each fit and adjust themselves according to the residuals at 
their respective energies. 
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Fig. 5. Combined 0.3-0.7 keV EPIC image of MCSNR 10512-6707. 
The red points delineate the contour level corresponding to 3(T above 
the average background surface brightness. The green solid line shows 
the best-fit ellipse to the contour, with the dashed lines indicating the 
Icr error on the fit. The white plus-sign marks the best-fit centre of the 
SNR. 


the resulting parameters do not provide a satisfactory physical 
interpretation of the X-ray emission. The addition of a second 
component in the form of a power law, provides the best fit to 
the spectrum, though again presents problems in the physical 
interpretation, which we now discuss. 


4.2.1. Soft emission 

We found no spectral signatures of ejecta, suggesting that the 
soft thermal plasma is dominated by swept-up ISM. This, and the 
low plasma temperature {kT ~ 0.2 keV), are consistent with an 
SNR in the Sedov phase of its evolution. Using the best-fit vapec 
parameters for both the single vapec and vapec-i-powerlaw 
models, and the Sedov solution, we estimated some properties 
of the remnant (see Sasaki et al. 2004; Kavanagh et al. 2015b, 
for example). However, for both the one and two component 
models, the derived SNR properties indicate a breakdown of 
the Sedov model with the determined upper limit on initial ex¬ 
plosion energy an order of magnitude less than the canonical 
10^* erg. Therefore, the soft thermal component of the MC¬ 
SNR J0512-6707 spectrum is not consistent with an SNR in the 
Sedov phase. 

If the remnant is not in the Sedov phase, then, given its rela¬ 
tively small size of ~ 25 x 22 pc, we might expect to see at least 
some ejecta contribution. One possible explanation that can ex¬ 
plain both the soft thermal emission of LMC abundance and the 
lack of ejecta emission is that the blast wave may have recently 
encountered a shell of dense material which it is now propagat¬ 
ing into. This would be the case if MCSNR J0512-6707 resulted 
from a CC explosion and the progenitor stellar wind has mod¬ 
ified the ambient medium creating a large, low-density cavity 
surrounded by a dense shell of ISM. Such stellar wind blown 
bubbles can have radii of the order of 10 pc (e.g., Weaver et al. 
1977; Chen et al. 2013), which is in agreement with the radius 
of MCSNR J0512-6707 of ~ 12 pc. In such a scenario, the rem¬ 
nant would have initially swept-up the circumstellar medium, 
then expanded into the low-density wind-blown cavity, before 
interacting with the dense shell (Dwarkadas 2005). When the 
blast wave encounters the dense shell, a shock is transmitted 
into the shell with significantly reduced velocity and a reflected 
shock travels back into the low-density interior which contains 
the ejecta. However, as shown by Dwarkadas & Dewey (2013) 
applied to the case of the Galactic remnant Kes 27, the emissiv- 
ity of the shocked ejecta is much lower than the shocked shell 
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Table 1. Spectral fit results MCSNR J0512-6707. See text for descrip¬ 
tion of the models. 


Component Parameter 

Value 

vapec 



phabs 

NH,Gal (10^2 Cm-2) 

0.05'' 

vphabs 

Nh,lmc (1022 cm^2) 

0.39 (0.32-0.48)'' 

vapec 

kT 

0.19(0.18-0.20)'' 


Z/Zo 

0.50 (fixed) 


EM (10^^ cm^2) 

1.0 (0.7-1.6) 


Fx‘^ (10^'“' erg s“' cm“2) 

2.9 


Lx"' (1024 ergs-') 

6.4 

Pit statistic 

Tv 

1.31 (736 d.o.f.) 

vapec -b powerlaw 

phabs 

NH,Gal (1022 cm-2) 

0.05" 

vphabs 

Nh,lmc (1022 cm-2) 

0.30 (0.23-0.41)'' 

vapec 

kT 

0.18 (0.15-0.19)'' 


Z/Zo 

0.50 (fixed) 


EM (lO^" cm-2) 

0.7 (0.5-1.5) 

powerlaw 

F 

2.07(1.81-2.26) 


norm (10-^) 

1.25 (1.01-1.59) 


Fx,totaf (10-'4 erg s-' cm-2) 

6.8 


Tx,vapec"'(1024 erg s-') 

3.9 


ix,powerlaw"'(1024ergs-') 

2.0 

Pit statistic 

Tv 

1.04(734 d.o.f.) 


Notes. The numbers in parentheses are the 90% confidence intervals. 

Fixed to the Galactic column density from the Dickey & Lockman 
(1990) HI maps. *** Absorption and thermal component abundances 
fixed to those of the LMC. 0.3-10 keV absorbed X-ray flux. 0.3- 
10 keV de-absorbed X-ray luminosity, adopting a distance of 50 kpc to 
the LMC. 


material, given the dramatic difference in densities. Therefore, 
only the shocked-shell emission would be observed which, given 
reduced transmitted shock velocity and the shell composition, 
would exhibit a soft thermal spectrum with ISM abundances. In 
Section 4.3, we refer to secondary evidence in the literature that 
MCSNR J0512-6707 does indeed result from a CC explosion 
and, therefore, the early evolution of the SNR into a wind-blown 
cavity is possible. 


4.2.2. Hard emission 

The origin of the hard emission presents a significant problem 
in its physical interpretation. We initially considered that the 
hard component may be due to a point source embedded in 
the extended emission of the SNR, most likely a background 
active galactic nucleus (AGN). During our analysis we per¬ 
formed a source detection using the SAS task edetect_chain 


(see Kavanagh et al. 2015b, for a full description of the de¬ 
tection method). However, the resulting sources located within 
the extent of SNR were identified as false detections due to the 
extended emission. We also searched in AGN catalogues for 
sources located within the extent of MCSNR JOS 12-6707 but 
found nothing. A further test of a possible AGN origin lies in the 
X-ray spectrum. If the hard component is the result of a back¬ 
ground AGN, we should expect that the hard component is more 
absorbed than the soft thermal component. Therefore, we per¬ 
formed spectral fits allowing for the additional absorption of the 
hard component. However, the hydrogen column density tended 
to zero, with an upper limit of ~ 1 x 10^' cm“^, indicating no ad¬ 
ditional absorption of the hard emission. For these reasons, we 
conclude that the hard X-ray emission is associated with MC¬ 
SNR JOS 12-6707 and not due to a point source. 

The hard emission is well fitted with either a simple 
power law or a high-kT, 1ow-t thermal plasma model of LMC 
abundance. Unfortunately, we are hampered in statistically 
distinguishing between these models because of the low count 
statistics at energies above ~ 2 keV. Therefore, we now discuss 
possible non-thermal and thermal origins of the emission. 

Non-thermal 

Morphologically, the SNR emission above ~ 2 keV is enhanced 
towards the northeast of the remnant, with the brightest hard X- 
ray emission correlated with the brightest region of radio emis¬ 
sion (see Fig. 6). In addition, this approximately correlates with 
the region where the polarisation fraction of the magnetic field is 
low and the magnetic field is oriented quasi-parallel to the SNR 
shock front (see Fig. 2), conditions favourable for particle ac¬ 
celeration and magnetic turbulence generation (Reynoso et al. 
2013). This suggests that the hard X-ray emission could be due 
to synchrotron emission from very high energy electrons. How¬ 
ever, shock speeds of the order of 10^ km s ' are still required 
for diffusive shock acceleration and X-ray synchrotron emission 
(see, e.g., Vink 2012). The temperature of the soft thermal X- 
rays is not compatible with such shock speeds. The reflected 
shock velocity may be high enough, especially relative to the ex¬ 
panding ejecta which it is moving into, though it is unclear if the 
synchrotron emissivity could explain the observed luminosity. 

To assess if the radio and X-ray data are compatible with 
the emission from a single electron population, we created a 
spectral energy distribution (SED). Integrated flux densities in 
the radio regime were reported in Reid et al. (2015, their Ta¬ 
ble 4). We adopt these values for the radio points on the SED. 
Eor the X-ray points, we used the EPIC-pn spectrum from 
Obs. ID 0741800201, the EPIC spectrum with the most counts. 
We then subtracted the background spectrum to ensure as much 
as possible that only X-rays due to MCSNR J0512-6707 were 
present. We confined our analysis to the 2-7 keV energy range 
to isolate the portion of the spectrum where the hard component 
is dominant. The radio and X-ray data points are shown in Pig. 7. 

We initially fitted the SED with a straight power law but 
it was immediately clear the X-ray fluxes were much lower 
than suggested by the radio data, which is expected if the 
electrons are close to the maximum energy of the relativistic 
electron distribution. We therefore applied a cut-off electron 
distribution of the form Ne(E) = where E is 

the electron energy, E^ax the cut-off energy, s is the spectral 
index of the electron distribution, and K a constant. The spectral 
index of the electron distribution is related to the spectral index 
of the photon distribution a as a = (s - l)/2. We make the 
simplifying assumption that each electron emits all its energy 
at its characteristic frequency (the (5-function approximation), 
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Fig. 6. XMM-Newton EPIC image of MCSNR 10512-6707 in the 2.0- 
4.5 keV range. The image was produced following the same procedure 
outlined in Section 3.2 and is overlaid with 20 cm radio contours from 
Fig. 3 - hottom left. 


and thus the resulting photon spectrum cuts off as e 
(Reynolds 1998). This cut-off function provides a very poor ht 
to the SED of MCSNR J0512-6707 with a = -0.61 (+0.06), 
and Vmax =1.8 (±0.8) x 10'® Hz (shown by the green dashed 
line in Fig. 7). Indeed, this fit is incompatible with the photon 
index confidence intervals determined in the hts to the hard 
X-ray emission (see Table 1). These confidence intervals are 
shown by the grey lines in Fig. 7. Therefore, we concluded that 
the combined radio and hard X-ray SED is not derived from the 
same electron population. This might be the case of the hard 
X-ray emission originates from specific regions of the SNR, as 
already indicated by its morphology. In the case of the SNR 
interaction with a dense shell, as suggested by the soft thermal 
X-rays, the integrated radio emission is made up of synchrotron 
emission from electrons behind the transmitted shock and any 
high velocity shocks, if present. The hard X-rays, on the other 
hand, are emitted near high velocity shocks only. We attempted 
to quantify the expected radio flux densities from electrons 
behind the high velocity shocks as follows: we adopted the 
radio spectral index of Reid et al. (2015) (a = -0.52), which 
is consistent with the typical value for SNRs (shown by the 
red dotted line in Fig. 7); we then hxed the value of a in 
the cut-off function to -0.52 and applied this model to the 
X-ray data only, resulting in the blue curve shown in Fig. 7. 
Comparing the blue curve to the radio data points indicates that 
the expected radio flux densities for a = -0.52 are more than 
two orders of magnitude lower than observed. However, we 
expect this to be a lower limit since young remnants with very 

fast shocks have radio spectral indices of a-0.7. Fixing a in 

the cut-off function to this value results in the magenta dash-dot 
line in Fig. 7, and the expected radio flux densities are now 
about one order of magnitude below the observed value, which 


is not an unreasonably low or high estimate. Unfortunately, 
without deeper and higher spatial resolution radio and X-ray 
observations, it is impossible to create SEDs for individual 
regions of the SNR and, therefore, to confirm if the hard X-rays 
result from synchrotron emission. 

Thermal 

The reflected shock may also be employed in a thermal origin for 
the hard X-rays. While the transmitted shock is slowed signih- 
cantly as it moves into the dense shell, the reflected shock travels 
back through the still expanding ejecta, further heating the al¬ 
ready shocked gas. Such a reflected shock-heating scenario has 
been used to explain the observed morphology at shock-cloud 
interaction regions in the Cygnus Foop (Fevenson et al. 2002). 

Another, though somewhat less likely, thermal explanation 
for the hard emission is the presence of a background galaxy 
cluster projected against the SNR. The hot intracluster medium 
of these objects can be described by CIE models with kT values 
up to ~ 10 keV (for a review, see Bohringer & Werner 2010). At 
such high temperatures, the spectrum is effectively featureless, 
which is especially true for low metallicity plasmas. Trial fits 
with a CIE plasma model in place of the power-law result in 
acceptable fits of equivalent goodness. However, the required 
chance alignment of a background galaxy cluster, located at the 
position of MCSNR J0512-6707 with a compatible morphol¬ 
ogy, makes this explanation for the hard component is extremely 
tenuous. In addition, the correlation of the synchrotron radio 
emission with Bremsstrahlung dominated galaxy cluster X-ray 
emission makes no physical sense. 

We have briefly discussed some possible origins for the hard 
X-ray emission from MCSNR J0512-6707 here but we cannot 
offer a dehnitive explanation. Deep high spatial resolution obser¬ 
vations with Chandra would allow the morphology of the hard 
emission region to be properly traced, helping to better identify 
the origin of the hard X-rays. 

4.3. Progenitor type 

The apparent association of MCSNR J0512-6707 with the 
DEM F97 HII region (Davies et al. 1976) and the proximity 
to the [WHO2011] A126 molecular cloud (Wong et al. 2011) 
would qualitatively favour a CC origin for the SNR. Unfortu¬ 
nately, given the lack of ejecta signatures in the X-ray spectrum, 
it is not possible to suggest a supernova remnant type directly. 

A secondary indicator for the SN type is the local stellar pop¬ 
ulation and star formation history (SFH). For example, a local 
stellar population harbouring many massive early-type stars and 
a corresponding recent burst of star formation would indicate 
that the SNR is most likely the result of a core-collapse explo¬ 
sion. Similarly, the absence of such a high mass stellar popula¬ 
tion and recent star formation burst would point to a Type la ori¬ 
gin. While this typing method is not infallible (see Maggi et al. 
2015, for a detailed description of the method and caveats) it 
still offers a good indication of the likely SN type and has been 
used in the past to gauge the explosion scenarios for a number of 
evolved EMC remnants (e.g., Kavanagh et al. 2013; Maggi et al. 
2014; Bozzetto et al. 2014). 

Maggi et al. (2015) performed a stellar population and 
SFH analysis of all confirmed SNRs in the EMC using pho¬ 
tometric data from Zaritsky et al. (2004) and the spatially re¬ 
solved SFH maps of Harris & Zaritsky (2009). They iden¬ 
tified 25 massive OB stars located within 100 pc of MC¬ 
SNR J0512-6707, inferring an expected CC/Type la ratio for 


Article number, page 9 of 11 





A&A proofs: manuscript no. astro-ph 



Fig. 7. SED of MCSNR J0512-6707. The data points at the lower fre¬ 
quencies are the radio flux densities reported in Reid et al. (2015). The 
data at the higher frequencies correspond to the background-subtracted 
EPIC-pn spectrum from Ohs ID 0741800201 in the 2-7 keV range, 
which excludes the contribution of the soft X-ray component. The red 
dotted line shows the power-law fit to the radio data, the green dashed 
line shows the exponentially cut-off power law fit to the entire SED, the 
magenta dash-dot line shows the exponentially cut-off power law fit to 
the X-ray data with a fixed to -0.7, and the blue solid line shows the ex¬ 
ponentially cut-off power law fit to the X-ray data with a fixed to -0.52 
determined by Reid et al. (2015). The grey lines indicate the upper and 
lower limit of the photon index in the X-ray power law fit (see Table 1). 


the region of 5.56^2 66 (their Table C.l). These factors indicate 
that MCSNR JOS 12-6707 most likely resulted from a CC event. 

5. Summary 

We presented a study of MCSNR JOS 12-6707 in the LMC and 
supplemented this with optical, radio, and IR data to gain an 
understanding of the physical properties of the SNR. The main 
findings of our analysis can be summarised as follows: 

- We showed that the ROSAT source [HP99 483] is the X-ray 
counterpart of the supernova remnant MCSNR JOS 12-6707 
identified by Reid et al. (201S). However, we determined 
that the remnant is larger than determined by Reid et al. 
(201S) from the optical emission. We suggested that the op¬ 
tical emission is confined to a region where the remnant ap¬ 
pears to be expanding into a denser ambient medium. 

- The 20 cm radio morphology exhibits a shell-like structure. 
Radio polarisation images at 3 cm and 6 cm in particu¬ 
lar, indicate a higher degree of polarisation in the northwest 
and southeast, tangentially oriented to the SNR shock front, 
indicative of an SNR where the expanding shell has com¬ 
pressed the magnetic field threading the ambient ISM. 

- The X-ray spectrum of MCSNR JOS 12-6707 is unusual in 
that it exhibits soft emission lines which are consistent with a 
soft thermal plasma in collisional ionisation equilibrium, but 
also shows featureless hard emission. We suggested that the 
soft thermal X-rays can be explained by a remnant initially 
evolving into the wind-blown cavity of a progenitor massive 
star and is now interacting with the surrounding dense shell. 


We discuss some possible thermal and non-thermal emission 
mechanisms for the hard X-ray component, though the true 
origin of these X-rays remains unclear. We found no spectral 
signatures of shocked ejecta. 

- The location of MCSNR JOS 12-6707 in the DEM L97 HII 
region and at the edge of the [WHO2011] A126 molecular 
cloud point to a core-collapse origin to the remnant. Due to 
the lack of ejecta signatures, we cannot directly determine 
the supernova type from the X-ray spectrum. Indirect evi¬ 
dence for the explosion mechanism is found in the study of 
the local stellar population and star formation history study 
of Maggi et al. (2015), which, because of the large number 
of nearby OB stars and recent star formation burst, suggest a 
core-collapse origin. This is consistent with the SNR evolu¬ 
tion into a the wind-blown cavity suggested above. 
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